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Abstract. We present Hubble Space Telescope Imaging Spectrograph long-slit spectroscopy of the emission line nebulae 
in the compact stee p spectrum radio sources 3C 67, 3C 277.1, and 3C 303.1. We derive BPT (Baldwin- Philips-Terlevich; 
iBaldwin etalJll98lh diagnostic emission line ratios for the nebulae which are consistent with a mix of shock excitation and 
photo-ionization in the extended gas. In addition, line ratios indicative of lower ionization gas are found to be associated with 
higher gas velocities. The results are consistent with a picture in which these galaxy scale radio sources interact with dense 
clouds in the interstellar medium of the host galaxies, shocking the clouds thereby ionizing and accelerating them. 

Key words, galaxies: active, galaxies: individual (3C 67), galaxies: individual (3C 277.1), galaxies: individual (3C 303.1), 
galaxies: quasars: emission lines. 



1. Introduction 

Powerful radio galaxies play a critical role in our understand- 
ing of both galaxy evolution and the phenomenon of activity in 
galactic nuclei. Yet we know little about how the radio galaxies 
are born and how they subsequently evolve. Recent work has 
identified the GHz Peaked Spectrum (GPS) and Compact Steep 
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Spectrum (CSS) radio sources as the most likely candidates 
for the progenitors of the large scale powerful classical double 
(FR2) sources (e.g.lO'Dea et al. 1 99 l llFanti et all 199(1 1 199.4 
iReadhead et alJll996blal lO'Dea & Baumlll997l) : for a review 
see lO'Deal Jl998h . The GPS and CSS sources are powerful 
but compact radio sources whose spectra are generally simple 
and convex with peaks near 1 GHz and 100 MHz respectively. 
The GPS sources are contained within the extent of the optical 
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narrow emission line region (< 1 kpc) while the CSS sources 
are contained within the host galaxy (< 15 kpc). 

Current models for the evolution of powerful radio galaxies 
suggest that these sources propagate from the ~ 10 pc to Mpc 
scales at roughly constant velocity through an ambient medium 
which declines in density as p(R) oc R 2 while the sources 
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densities of powerful radio sources as a function of linea r size 
(from tens of parsecs to hundreds of kpc) (e.g. lO'Dea & B aum 
ll997tlFanti etaTfcoOll) . However, the situation must be more 
complicated than this simple picture. We give two reasons: 
(1) The GPS and CSS sources must interact with the host 
galaxy as they propagate through it. The discovery of emission 
line gas aligned with and presumably co-spatial with the 
CSS radio sources and the pres ence of broad and co mplex 
integrated emission line profiles jGelderman & Whittldfl99l 
indicates that the radio source is st r ongly interacting with th e 
ambient gas Jde Vries etaP 1 19971 1 19991: lAxon et alJ l2000l) . 
Therefore, we would expect s hocks to contribute s trongly to 
the ionization of the gas (e.g. iBicknell et alJll997l) . (2) The 
GPS sources are observed to have expansion velocities several 
times higher (Conway 1998, private communication) than 
the estimated advance spee ds of large scale classical doubles 
^Alexander & Leahvlll987l) . This would require the evolving 
GPS sources to decelerate as they propagate though the host 
galaxy and would require the radio sources to dim faster than 
the simple models predict. It may be that the deceleration takes 
place via interaction with ambient gas (see iDe Younglll993t 
IC a rvalholl994fl99l). 

We are carrying out a study of the kinematics and ion- 
ization of the aligned emission line nebulae in three CSS 
radio sources: 3C 67 (galaxy, z = 0.310, linear size D=10.1 
Kpc), 3C 277.1 (quasar, z = 0.321, D=6.9 Kpc), 3C 303.1 
(galaxy, z = 0. 267, D=6.2 Kpc) 1 with HST/STIS long-slit 
spectroscopy. In lO'Dea et alJ J2002I) we presented our results 
on the kinematics of the [O m] /15007 emission line and found 
complex emission line profiles and large differences in velocity 
offset on the two sides of the nucleus, suggesting that the 
cloud motions are being driven by shocks i nduced by the 
expanding radio lobes. In lO'Dea et alJ (|2003) we discussed 
two models for the cloud kinematics - (1) acceleration by the 
radio source bow shock and (2) acceleration by the post bow 
shock wind, and concluded that the bow shock acceleration 
was favored. Here we present the results of our low dispersion 
spectra of several diagnostic emission lines. We compare our 
results with shock and photo-ionization models and examine 
the relationships between the ionization diagnostics and the 
gas kinematics. 



1 We adopt a Hubble constant of H„ = 75 km s 1 Mpc 1 and a 
deceleration parameter of q„ = 0.0. 



2. Observations 

We obtained STIS long slit spectra through the 52" x 0.1" 
slit through several gratings. The instrumental parameters are 
summarized in Tabled The observations through the medium 
dispersion G750M grism has been used primarily to study the 
kinem atics of the nebula and are discussed by lO'Dea et alJ 
d2002h . For the study of the ionization, we observed through 
the low dispersion G750L (5236 to 10266 A) and G430L (2900 
to 5700 A) grisms centered on 7751 A and 4300A, respectively, 
with a spectral resolution of 5A/pixel and a Line Spread 
Function (LSF) of 2.0 pixels resulting in a velocity resolution 
of ~ 700 km/s at the center of the blue side of the spectra and 
~ 390 km/s at the center of the red side 2 . The spatial resolution 
is given by the Point Spread Function (PSF) of the detector, 
with a FWHM of 2.3 pixels. 

We integrated for one orbit on each slit position (~ 
2500-3000 sec) with the orientation taken roughly parallel 
to the radio source axis. In each source one slit was centered 
on the source nucleus and aligned along the radio axis. In 
3C 67 and 3C 303.1, where the emission line gas is slightly 
misaligned with respect to the radio axis, we rotated the slit 
to place it along the emission-line nebula axis, rather than the 
radio source axis. 



3. Data reduction 

The standard STIS reduction pipeline was used to remove 
detector signatures such as bias, dark current, and flat fielding 
and to apply the flux calibration. Cosmic ray hits can be 
quite numerous over the course of orbit-long exposures. We 
therefore split each exposure per orbit into two equal-length 
parts in order to allow removal of cosmic rays. Any surviving 
cosmic rays were removed and replaced by the average of the 
flux values in the adjacent pixels around the cosmic ray. 

To correct the spatial dispersion of the light produced 
by the spectrograph (given by the PSF) we have averaged 
(weighted by the flux errors given by the spectrograph) 
three pixels in the spatial direction for every source. We also 
averaged all the rows in each lobe for each source, to produce 
an average lobe spectrum for each side of the source. 



We used the IRAF/STSDAS Specfit dKrisslll994 software 
to fit Gaussians to the profile of each emission line in our 
sources (Tableland Fi gures 17721 and 01 measuring each line's 
integrated flux, full width at half maximum (FWHM) and 
central wavelength (i.e., velocity offset relative to the nucleus). 
Some constraints were adopted to reduce the number of free 
parameters and limit our options to "physically consistent 
models", as follows: We classified the detected lines in the 
spectrum as either high or low ionization (see Table [2j. We 
required all gas with same ionization state to have the same 
FWHM and velocity offset, and took as as free parameters the 
widths and velocities of H/3 and [Oin]/15007. The integrated 

2 Corrected for redshift. 
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fluxes of every line were free to vary exc ept for those with 
known ratios (from atomic physics, see e.g.. lOsterbrockll9 89). 
in this case: the [Nn] XX 6548,84 and [Om] XX 4959,5007 
doublets. We used a power law to fit the continuum. Whenever 
this model was not accurate enough in parts of the spectrum, 
we limited the fit to the adjacent regions of each emission 
line, improving our continuum model and producing accurate 
fits for the line profiles. Only a few of the potential emis- 
sion lines (Table [5J were detected in the extended nebulae 
(signal to noise ratio > 3) - the [On]/i/l 3727,29 doublet, 
Ha, U/3, [Om] XX 4959,5007 and the [Nn] XX 6548,84 , 
[Sn] XX 6716,31 doublets 3 . Upper limits were defined to 
be three times the noise of the spectrum, multiplied by 
the width of the line: 3xRMSxFWHM. We used the width 
we measured in the resolved lines of the same ionization status. 

Specfit is known to overestimate the errors for flux and 
central wavelength measurements in our data. To improve 
the estimation of our error bars we created over one hundred 
artificial spectra consisting of a perfect Gaussian with known 
random white Gaussian noise added, which covered the 
parameter space of our data in signal to noise ratios and 
resolution. We created ten artificial spectra for each different 
signal to noise ratio and ran Specfit on each spectrum. We 
compared the Specfit derived errors with those expected based 
on the parameters of the artificial spectra and estimated a 
correction to the Specfit errors. The final step in the data 
reduction was correcting the data for reddening, u sing the 
Galactic dereddening curve in ICardelli et alJ \ 19891) and t he 
measured Galactic extinction values of lSchlegel et alJ dl998h . 

lO'Dea et alJ (2002) searched for the possible presence 
of faint, broad wings (> 1500 - 2000 km/s) using the low 
dispersion data, that might have been undetected in the higher 
dispersion data. They found no strong evidence for broad, 
non-Gaussian components substantially above ~ 1500 km/s in 
any of these galaxies. 



4. Results 

We have detected several lines in the extended nebulae of 
all three CSS sources - [On], [Om], H/3, Ha, and [Nn]. We 
have used these lines to determine the standard diagnostics 
of ionization and compared them to the predictions of shock 
and photo-ionization models. We have also compared the 
ionization le vels with the gas kinematics, as presented and 
discussed bv lO'Dea et al] J2002I) . 



4.1. Ionization diagnostics 

Given our line detections we were able to construct one 
of the key diagnostic Baldwin-Philips-Terlevich (BPT) di- 

3 We have detected the [S n] AA 6716,31 emission in very few points 
of the sources and our low wavelength resolution does not allow us to 
deblend the doublet components. 



agrams Jialdwin et alJ 1 198 it IVeilleux & Osterbrockl Il987t 
Mov & Roc^a-Volmeraiigell2002l) . namely the fOnl/fOml vs. 
[Om]/H/? diagram (Figs. l4"l5lol Other potentially useful BPT 
diagrams (e.g., [Sn]/Ha vs. [Nn] 6583/Ha, [On]XX 3727+29/ 
[Om]5007vs. [Nn] 6583/Ha, [Om] 5007/H/3 vs. [Sn]/Ha, 
etc.) were not possible given the limited number of detected 
lines. We have compared our results with the two main classes 
of ionizat ion mechanisms: shock i onization usin g MAPPINGS 
III ( JKewlev et alJ 12003 in preol iDooita & Sutherland 1 19961) 
and AGN photo-ionization using CLOUDY ( Ferland et al. 
Jl998t) \ The MAPPINGS shock code calculates the spectra 
from gas which has been directly ionized by the shock (pure 
shock models) and unshocked gas which has been ionized by 
the radiation emitted by the hot shocked gas (pure precursor 
models). In our BPT diagrams (Figs. I4I5I6> we show the 
predicted line ratios calculated by MAPPINGS for both a pure 
precursor and pure shock model, as well as linear combinations 
of these two (30%, 50%, 70% contribution to the emission 
line luminosity from shocks), including a range of magnetic 
fields (MAPPINGS 's magnetic parameter from 0.5 to 10) 
and shock velocities from 100 to 1000 km/s. In addition we 
include the predicted line ratios calculated by CLOUDY for 
a range of ionization parameter (Log U ranging from -2 to 
-3.8) and cloud density (Log n = 0,2,4). We note that "real" 
radio galaxies are probably more complicated than assumed in 
the MAPPINGS and CLOUDY calculations and suggest that 
the results be taken with caution. 

We find that the data points tend to spread across the 
diagrams (Figs. I4I5I6> covering about one dex in both di- 
mensions. The nuclei lie to the lower right of the distribution 
of points (with the highest values of [Om]/H/3, and lowest 
values of [On]/[Om]) consistent with the nuclei containing 
the highest ionization gas. The data for the extended emission 
clearly indicate lower ionization than the nuclei. In addition 
the data for the extended emission are at higher ionization 
than the 100% shock models from MAPPINGS and are at 
lower ionization than the CLOUDY AGN photo-ionization 
models. In general the extended emission lines are consistent 
with a mixture of shocked and photo-ionized gas. In 3C 67, 
the data lie between the MAPPINGS models for a contribution 
to the observed luminosity from shocked gas ranging from 
to 50%. In 3C 303.1, the contribution to the luminosity 
from shocked gas is between 30 and 70%. In 3C 277.1, 
the data scatter around the model for 100% contribution to 
the luminosity from precursor gas. Thus, the contribution 
from shocks increases from 3C 277.1 to 3C 67 to 3C 303.1. 
The sources tend to lie in the regions for moderate to high 
shock velocities (500 to +1000 km/s) and we do not obtain 
any useful constraints on magnetic field strengths in the 
MAPPINGS models. We note that these results are consistent 
with previous work on CSS radio sources using ground b ased 
data JGelderman & Whittle! fl994l iMorganti et all 1 19971) . We 
also used the observed Ha/H/3 ratios to test the hypothesis that 
the low values of the [O n]/[0 m] ratio are produced by large 
amounts of redenning. However, the observed Balmer ratios 
(or even values 10% higher than observed) are not consistent 
with significant redenning effects on the [Oii]/[Om] ratio 
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(except possibly in the nucleus of 3C277.1). 



4.2. Can the central AGN power the emission line 
luminosity in the extended nebulae? 

We compared the number of ionizing photons produced by 
the nucleus of the source, with the number of photons needed 
to produce the observed emission line luminosity (see e.g . 
Wilson et alll988|lBaum & Heckmanll989llAxon et all2 000: 
O'D ea et al J 12000). Assuming radiative recombination under 
case B conditions, the number of ionizing photons Afyj needed 
to produce the observed H/? luminosity is: 

Nufj = 2.1 x lO 12 ^ photons s" 1 

We use the integra te d rO ml/t5007 fluxes measured by 
iGelderman & Whittle! ( 1 1994 (see table and scale us- 
ing the typical ratio for t he narrow line component s 
in CSS sources (see also IGelderman & Whittle 1994): 
Hf3/[Om] ,15007 = 0.18 + 0.02. 

The number of photons in the continuum, between fre- 
quencies v\ and V2 are given by: 

AW = 4nD 2 S (ahy l (v° x - vf ) photons s _1 

where D is the luminosity distance, the flux density spectrum 
is given by F v = Sov~ a (we adopt a-l) and h is Planck's 
constant. We are only interested in the photons with enough 
energy to ionize Hydrogen, so we choose v\ = 3.3xl0 15 Hz 
(912A or 13.6eV) and v 2 = 4.8xl0 17 Hz (17 keV). For our 
spectral index, a=\, higher frequencies do not add a significant 
number of photons. Note that this analysis is subject to 
the caveat that the continuum emission may not be emitted 
isotropically, and the exte nded nebulae may see a different 
luminosity than we do (e.g. lPenston et alJl990t) 

The results are shown in Table |4] We find that the nu- 
cleus apparently produces enough ionizing photons to power 
the emission line luminosity in 3C 277.1 and possibly 3C 67, 
but not 3C 303.1. This is consistent with the results from our 
BPT diagrams. 

4.3. Diagnostic ratios vs. kinematics and distance 

Here we look for additional clues to the nature of the ionization 
and acceleration of the emission line gas. In shock models, 
the ionization sh ould be related to the gas kinematics (e.g., 
IClarketalJl998l while if the ionization is dominated by AGN 
photoionization, the ionization level and kinematics should be 
largely independent. 

We have studied the behavior of the diagnostic ratios 
[Om] 5007/H/J, [On] ^3727+29 / [Om] ^5007, [Nn] 6583/ 
Ha, [S n]/Ha; with distance to the nucleus of the source, and 
the gas kinematics (velocity offset relative to the nucleus, 



and FWHM). The kinematic data have been ta ken from the 
medium dispersion results rtO'Dea et alJ 120021) given their 
superior spectral resolution. The low dispersion data give 
consistent results. The relationships are plotted in Figures Q to 
[8] and the results are summarized in Table [5] We have applied 
a Spearman rank test to search for correlations and regard 
correlations with a significance >99% as real. 

The three sources show a decrease in ionization with 
distance from the nucleus (based on [Oiii] 5007/HjS, 
[On]/U 3727+29/[Oiii]5007). This behavior (high ion- 
ization in the nucleus, low in the lobes) can be explained by (1) 
the influence of the central engine in the inner regions, and a 
low ionization source, like shocks, further out, and (2) dilution 
of photons from the nucleus with distance. The sources 
show a weak trend for ionization to decrease with increasing 
FWHM and a strong trend for ionization to decrease with 
increasing velocity offset. This relation between ionization 
level and the kinematics is consistent wi th shock ionizatio n 
and acceleration of the gas (as suggested bv lO'Dea et alJ2002l) . 



4.4. Cloud properties 

We have estimated the electron temperature using the 
[Om]/14363/[Ora]/l5007 diagnostic lines. Given our limited 
signal-to-noise we derive an estimate of the <?~ temperature in 
the nucleus of 3C277.1 of ~20,000K and obtain lower limits 
(>7,000K) for the rest of the sources. These temperatures 
are consistent with previous measurements of temperatures 
in the nucleus and extended emission-line regions of A GN 
dClark et alJll998HWilson et alJll997tlStorchi-Bergmann et alJ 
1996). Tadhu nter et alJ (119891) showed that the observed ratio 
O iii/14363/O ni/15007 implied for several extended emission- 
line regions electron temperatures ranging from 12800K to 
22000K, whereas photoionization models predict much lower 
values. They discussed several possible explanations for this 
discrepancy between the results and the models, including 
additional heating sources such as cosmic rays or shocks, or 
metal abundances lower than the solar values assumed in the 
models. 

In 3C 277.1 we also found high Ha/H/Jvalues for the 
nucleus (10.4 + 0.2) consistent wi t h the measurements of 
CSS quasars by iBaker & Hunsteadi dl995l) . who argued for 
very dusty gas in these objects. On the other hand, 3C 67 and 
3C 303.1 show Ha/Hj8 ratios of ~ 3.1. 



5. Summary 

We present Hubble Space Telescope Imaging Spectrograph 
long-slit spectroscopy of the aligned emission line nebulae in 
three compact steep spectrum radio sources: 3C 67, 3C 277.1, 
and 3C 303.1. In previous papers we have reported evidence 
that the kinematics of the gas is consistent with being driven 
by shocks from the expanding radio lobes. Here we present 
the results of our low dispersion spectra of several diagnostic 
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emission lines. We compare our results with shock and photo- 
ionization models and examine the relationships between the 
ionization diagnostics and the gas kinematics. 

We find that the ionization diagnostics are consistent with 
a mix of shock and photo-ionization in the extended emission 
line gas. The data for the extended emission are at higher ion- 
ization than the 100% shock models from MAPPINGS and are 
at lower ionization than the CLOUDY AGN photo-ionization 
models. In general the extended emission lines are consistent 
with a mixture of shocked and photo-ionized gas. In 3C 67, 
the data lie between the MAPPINGS models for a contribution 
to the observed luminosity from shocked gas ranging from 
to 50%. In 3C 303.1, the contribution to the luminosity 
from shocked gas is between 30 and 70%. In 3C 277.1, the 
data scatter around the model for 100% contribution to the 
luminosity from precursor gas. The sources tend to lie in the 
regions for moderate to high shock velocities (500 to 1000 
km/s) and we do not obtain any useful constraints on magnetic 
field strengths in the MAPPINGS models. 

The three sources show a decrease in ionization with 
distance from the nucleus (consistent with a decrease in 
photoionization with distance) a weak trend for ionization 
to decrease with increasing FWHM and a strong trend for 
ionization to decrease with increasing velocity offset (which is 
consistent with shock ionization) . 

These results are consistent with a picture in which the CSS 
sources interact with dense clouds as they propagate through 
their host galaxies, shocking the clouds thereby ion izing and 
accelerating them (as suggested bv lO'Dea et al.l2002l) . 
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Fig. 1. Spectra of the lobes and nucleus of 3C 67. Left is G430L 
and right is G750L. The top row corresponds to the spatially 
averaged spectra of the southern lobe, the middle row to the 
averaged spectra of the northern lobe and the bottom row cor- 
responds to the nucleus. We have trimmed the sides of some 
spectra for clarity but all the wavelengths of interest are shown. 
The units of flux are erg s _1 cm _2 arcsec _2 A -2 . 



Table 1. Journal of HST/STIS Observations. 



Source 


Grism 


Central A 


PA 


Time 






A 


Degrees 


Sec 


3C 67 


G750L 


7751 


17 


2160 


3C 67 


G430L 


4300 


17 


2660 


3C 277.1 


G750L 


7751 


311 


2340 


3C 277.1 


G430L 


4300 


311 


2800 


3C 303.1 


G750L 


7781 


331 


2900 


3C 303.1 


G430L 


4300 


331 


5300 




Wavelengths (angstroms) Wavdcrisjrti i: t L 'jrr,Tn-i 



Fig. 2. Spectra of the lobes and nucleus of 3C 277.1. Left is 
G430L and right is G750L. The top row corresponds to the 
spatially averaged spectra of the southern lobe, the middle row 
to the averaged spectra of the northern lobe and the bottom row 
corresponds to the nucleus. We have trimmed the sides of some 
spectra for clarity but all the wavelengths of interest are shown. 
The units of flux are erg s _1 cm _2 arcsec _2 A -2 . 

Table 4. Photon counting. 



Source 


Distance 


Log N Hf! 


Log N w „ c 


N H /j/N Nuc 


3C 67 


1620 


54.4 


54.3 


1.3 


3C 277.1 


1700 


54.6 


55.1 


0.35 


3C 303.1 


1350 


54.4 


53.1 


17.0 



Source, distance to the source in Mpc (Ho=75km/s/Mpc) and 
number of photons needed to ionize H^, ionizing photons produced 
by the nucleus, and the ratio between the last two. If the ratio is <1, 
the nucleus is producing enough photons to ionize tig, if the ratio 
is higher than one, another source of ionization, such as shocks, is 
required. 



HST proposal ID=8104. The slit was placed parallel to the 
radio source axis. 
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Fig. 3. Spectra of the lobes and nucleus of 3C 303.1. Left is 
G430L and right is G750L. The top row corresponds to the 
spatially averaged spectra of the southern lobe, the middle row 
to the averaged spectra of the northern lobe and the bottom row 
corresponds to the nucleus. We have trimmed the sides of some 
spectra for clarity but all the wavelengths of interest are shown. 
The units of flux are erg s _1 cm~ 2 arcsec _2 A 2 . 




Log ( [01 1 1] 5007 / H&N) 



Fig. 4. 3C 67. Reddening-corrected 

[On] Xt3727+29/[Oiii] ,15007 vs. [Om] ^5007/H^ inten- 
sity ratios. Symbols for the models are shown in key to the 
figure. Data is shown for the nucleus (Nuc) and at each 0.15 
arcsec along the slit in the extended region, labeled with their 
distance to the nucleus. We also show the points for the aver- 
aged extended emission in the northern (N) and southern (S) 
lobes. The arrows in the plot represent those points where lines 
were not detected, using an upper limit of (3*RMS*FWHM). 
We compare our results with the t wo main classes of ionizatio n 
mechanisms: MAPPINGS III f llKewlev et alJ 12003 in prepl : 
Dopita & Sutherland 1996) shock models (solid lines) and 
CLOUDY dFerland et alJ i 19981) ) AGN photo-ionization mod- 
els (dashed lines). For the shock models from MAPPINGS 
we included both pure precursor and pure shock models, as 
well as linear combinations of these two (30%, 50%, 70% 
contribution to the emission line luminosity from shocks). The 
results do not depend on magnetic field strength, and so for 
simplicity we show models with B = 10. Shock velocity in 
this models ranges from 100 to 1000 km/s and increases to 
the right. For the CLOUDY AGN photoionization models we 
include a range of ionization parameter (Log U ranging from 
-1.8 to -3.8) and cloud density (Log n = 0,2,4). 




Lag ([QUI] 5007/ HftN) 



Fig. 5. Same as figuregl but for 3C 277.1. 
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Fig. 6. Same as figure^ but for 3C 303.1. 
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Fig. 7. Behavior of the diagnostic ratio [On]/U 3727+29/ 
[Ora]5007versus distance, velocity offset and FWHM. 
Diamonds represent 3C 67, squares 3C 277.1 and triangles 
3C 303.1. 
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Fig. 8. Behavior of the diagnostic ratio [Om] 5007/H/3 ver- 
sus distance, velocity offset and full width half maximum. 
Diamonds represent 3C 67, squares 3C 277.1 and triangles 
3C 303.1. 
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A (A) 






3C 67 






3C 277.1 






3C 303.1 




Line 


Ion. 


Nuc 


S 


N 


Nuc 


S 


N 


Nuc 


S 


N 


[Mgn] 


2799.12 


Low 


24.0±0.7 


- 


- 


2738±40 


- 


- 


- 


- 


- 


[Nev] 


3425.90 


High 


143.9±6.3 


23.2±1.4 


- 


427±18 


- 


- 


- 


- 


- 


[On] 


3727.37 


Low 


111.2±4.6 


57.1+9.8 


20.3±1.5 


188.9+15 


65.0±1.8 


33.8±0.6 


275.4±3.6 


144.2±1.4 


63.1 + 1.3 


[Nera] 


3868.76 


High 


133.6±6.0 


- 


- 


627±17 


- 


- 


52.4±1.8 


- 


- 


[Nem] 


3967.47 


High 


13.2+4.1 


- 


- 


412±26 


- 


- 


- 


- 


- 


HS 


4101.73 


Low 


- 


- 


- 


212.8±3.3 


- 


- 


- 


- 


- 


Hy 


4340.47 


Low 


- 


- 


- 


343.1±4.0 


- 


- 


25.1+1.1 


- 


- 


[Om] 


4363.21 


High 




15.4+1.1 




332.7±3.4 












[He ii] 


4685.68 


Low 








57.7±3.7 






23.7+1.1 






w 


4861.33 


Low 


95.9±2.7 


33.4+1.4 




281.0±4.1 


30.4±1.6 




72.7±1.5 


13.7±0.6 


18.9±0.8 


[Om] 


4958.92 


High 


54.0±3.2 


67.9±1.2 


25.2±0.9 


1031.3+4.7 


27.4±1.2 


35.7±0.8 


258.9±1.6 


32.3±0.6 


35.9±0.6 


[Om] 


5006.85 


High 


1566.8±5.7 


196.7±1.4 


73.0±1.0 


2998.3±9.2 


79.7+1.3 


103.9±0.9 


752.1±2.8 


95.3±0.6 


104.4±0.7 


[He i] 


5875.70 


Low 


45.2+3.1 






98.4±4.5 












[Oi] 


6300.31 


Low 




















[Nn] 


6548.06 


Low 








181.6+5.1 






59.6±2.7 


16.8+0.8 




Ho- 


6562.82 


Low 


372.0±3.7 


103.4+1.7 


37.6+1.3 


2923.3±12.6 


94.1+2.5 


88.7±2.0 


210.6±2.3 


39.7±0.9 


59.4+1.1 


[Nn] 


6583.39 


Low 


75.2±3.3 


58.9+1.6 




527.4±5.9 


4.9±2.4 


36.1±2.0 


193.5+3.1 


48.9±0.9 


34.5±1.0 


[Sn] 


6716.42 


Low 








160.6±5.2 






56.2+16 


32.7+1.1 




[Sn] 


6731.78 


Low 





















List of the emission lines we searched for in our sample. The first column gives the s pecies; the second column the rest-frame wave- 
length; the third column the ionization status. For wavelength and ionization references see Eracleous & Haloern 1 2004). For each source we 
present the integrated flux (in 10~ 16 erg s~'cirr 2 arcsec~ 2 ) for the detected lines in the nucleus (Nuc) and southern (S) or northern (N) lobes. 
Nondetections are indicated by For the lobes we have used the spatially averaged spectra. However, a line marked in the table as a non 
detection in the averaged spectra might have been detected in a couple of very localized positions. We require 3<r for a detection. See also 
FiguresQ[3]and[3]for spectra of the sources. 



Table 3. Source Properties 



Parameter 


3C 67 


3C 277.1 


3C 303.1 


ID 


G 


Q 


G 


redshift 


0.310 


0.321 


0.267 


scale (kpc/arcsec) 


4.04 


4.13 


3.65 


radio power logioPso/y; (Watts Hz -1 ) 


26.3 


26.4 


26.0 


angular size 6 (arcsec) 


2.5 


1.67 


1.7 


linear size D (kpc) 


10.1 


6.9 


6.2 


integrated emission line flux F(Ora /15007) (10~ 15 ergs s _I cirT 2 ) 


22 


31 


28 


integrated ([Oni/15007) line width at 50% of line peak (km s~') 


600 ± 42 


510 ± 19 


815 ± 18 


integrated (Om /15007) line width at 20% of line peak (km s~') 


915 ±85 


740 + 50 


1170 ±41 


spectral age (yr) 


5 x 10 4 


2 x 10 s 


1 x 10 5 


advance speed (v/c) 


0.22 


0.04 


0.07 



We adopt a Hubble consta nt of H = 75 km s 1 Mpc 1 and a deceleration parameter of q a = 0.0. The integrated emission line flux 
and width are from lGelderman & Whittle] J 1994 . The spectral age is estimated by fitting a continuous injection model to the integrated radio 
spectrum and is taken from lMurgia^HrtT P^^ The advance speed is estimated using 2 v = linear size / spectral age. 
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Table 5. Behavior of diagnostic ratios with distance and kine- 
matics. 



Source 


Plot 


[Om] 5007/H/3 


[O II]/[0 III]" 


[Nil] 6583/Ha 


[S n]/Ha 




Dist 


Decrease (99%) 


Increase (95%) 


Asymmetric^) (95%) 


No trend 1 


3C 67 


VO 


Decrease 2 (99%) 


Increase 2 (85%) 


Increase 2 (95%) 


No trend 1,2 




FWHM 


No Trend 


No Trend 


No Trend 


No trend 1 




Dist 


Decrease 4 (95%) 


Increase (99.5% S) 


Decrease (90% N) 5 


No trend 1 


3C 277.1 


VO 


Decrease 4 (90%) 


Increase (90% S) 


Increase(95% N) 5 


No trend 1 




FWHM 


No trend 4 


Increase (95% S) 


No trend 


No trend 1 




Dist 


Decrease (99.5%) 


Increase (99%) 


No trend 3 


Asymmetric 3 (99.5%) 


3C 303.1 


VO 


Decrease (99%) 


Increase (99.5%) 


No trend 


Increase (95%) 




FWHM 


Increase (90%) 


No Trend 


No trend 


No trend 



<a)[0 II],U3727+29/[0 III]i5007 

•Very few clear detections 

2 Most of the points around \ VO\ -150 km/s 

3 Increases to the south. Decreases to the north. 

4 Few clear detections of BfS 

5 Few clear detections of [N n] in the South Lobe. 

Summary of the behavior of the different ratios with distance to the nucleus, velocity offset (relative to the nucleus) and FWHM. First 
column gives the name of the source, second column the "X" axis in the plot; the next four columns give the studied ratios. Asymmetric: 
Significantly different values in the two lobes. No trend: No observable trend. Increase/Decrease: Trend of the ratio with increasing values 
of distance, velocity offset or FWHM. Where trends were found, we note the -minimum- confidence level based on the Spearman's Rank 
Correlation for that plot. A letter S or N means the correlation was only found in the southern (S) or northern (N) lobe. See the text for more 
detailed comments on each ratio and source. 



